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Summary 

The objective of this study is to predict the subsonic 
flow field around several High Speed Research Refer- 
ence H configurations at various pitch and yaw angles. 
The configurations achieve a high lift-to-drag ratio with 
leading- and trailing-edge flap deflections that are opti- 
mized to minimize flow separation around the leading 
edge and flap hinge lines. A sequence of structured- 
viscous grids was generated for the high-lift Reference H 
configuration. The first grid modeled the wing-body 
high-lift geometry, and the second grid incorporated the 
nacelles and the horizontal tail. The third grid modeled 
the full-span geometry for sideslip calculations, and was 
obtained by mirroring a coarser version of the second 
grid. The CFL3D code, a Reynolds averaged, thin-layer 
Navier-Stokes flow solver for structural grids, was used 
for the flow solver and modeled the free-air Reference H 
high-lift configuration at wind tunnel conditions of Mach 
number 0.24 and Reynolds number of 1.4 x 10 5 per in. 
Pitch sweeps were performed at angles of attack from 6° 
to 15°. Sideslip angle sweeps at 0° < |3 < +18° were per- 
formed at an angle of attack of 8°. 

The lateral and longitudinal performance character- 
istics were well predicted and very good force and 
moment comparisons are seen. The pressure distributions 
show generally good comparison and trend at the lower 
angles of attack; however, comparisons diverge when the 
predicted flow on the outboard portion of the wing 
becomes extensively separated, which occurs at about 
12° angle of attack. A very complex multiple vortical 
system is predicted at the higher angles of attack. Also, at 
the higher sideslip angles, an extensive area of stalled 
flow forms on the left side of the vertical tail. Postpro- 
cessing of the numerical solutions provided detailed sur- 
face and off-surface flow quantities, not obtained in the 
wind tunnel. The plotted off-surface pressure contours of 
the computational vortical flow field was very insightful 
in interpreting the wind tunnel data and provided a com- 
prehensive three-dimensional understanding of the flow. 

Introduction 

The development of technologies for an economi- 
cally feasible high speed civil transport by NASA and the 
U.S. aerospace industry is currently underway, and is the 
main focus for the High Speed Research program. Simi- 
lar to the Supersonic Cruise Research program, which 
was initiated in July 1972 by NASA, the High Speed 
Research program has focused on the potential benefits 
and trade-offs of advancements in aerodynamic effi- 
ciency, structures and materials, propulsion systems, and 
stability and control requirements applied to advanced 
supersonic cruise aircraft concepts. From a performance 
point of view, the aircraft would be designed for high 


aerodynamic efficiency in supersonic cruise and have 
good low-speed characteristics. For efficient supersonic 
cruise, the aircraft must have highly swept, low-aspect- 
ratio wings; however, this type of configuration is ineffi- 
cient in the subsonic flight regimes. Dining take-off and 
approach, the flow separates from the leading edge of the 
wing and forms a vortex system, which causes an 
increase in drag. The additional thrust required from the 
engines to overcome the increased drag and the low 
speed performance deficiencies, generates an unaccept- 
able high level of engine noise. The environmental 
impact of the engine-generated noise can be minimized 
by reducing flow separation, which maximizes the take- 
off lift-to-drag ratio. The environmental impact is mini- 
mized by incorporating special leading- and trailing-edge 
devices or attached-flow flaps onto the aircraft. 

The focus of this investigation is the accurate 
prediction of the subsonic flow field about a Reference H 
model with leading- and trailing-edge flaps using a 
proven viscous numerical method. Computations were 
performed at various pitch and yaw angles using 
structured viscous grids for the wing-body, for the wing- 
body with nacelles and horizontal tail, and for the wing- 
body with nacelles and horizontal and vertical tail 
configurations. 

Symbols and Abbreviations 

b wing span, in. 

Q> drag coefficient, Drag/<?„5 

C L lift coefficient, HixJq^S 

Ci rolling-moment coefficient (positive right 

wing down). Rolling momentA^Sb), body 
axis 

C m pitching-moment coefficient about 0.50c 

(positive nose up), Pitching momentA^^c) 

pitching-moment coefficient at zero lift 

C„ yawing-moment coefficient (positive nose 

right). Yawing moment/(^ 00 5&), body axis 

C p static pressure coefficient, {p - p M )/q x 

Cy side force coefficient (positive right), 

Side totcd(q^S) 

c mean aerodynamic chord, in. 

q local chord, in. 

ij,k grid index 

M ffee-stream Mach number 

p static pressure, psi 

p t total pressure, psi 

/?„ ffee-stream static pressure, psi 

q m free-stream dynamic pressure, psi 



Re^ Reynolds number per unit length, L/in. 

S reference area, in 2 

V„ free-stream velocity, in/sec 

xi local x coordinate, zero at leading edge, in. 

x/C[ fraction of local chord 

x,y,z body axis coordinates, in. 

y + law-of-the-wall coordinate 

y/(b/2)i fraction of local semispan 
a angle of attack, deg 

P angle of sideslip, deg 

6le leading-edge flap deflection angle (positive 

down), deg 

§xe trailing-edge flap deflection angle (positive 

down), deg 

curvilinear coordinates 
Abbreviations: 

configl wing-body configuration 

config2 wing-body with nacelle and horizontal tail 
configuration 

config3 full-span wing-body with nacelle and horizon- 
tal and vertical tail configuration 

LE leading edge 

TE trailing edge 

MBL model buttline station, in. 

MFS model fuselage station, in. 

MW megaword 

Test and Model 

The NASA Langley 14- by 22-Foot Subsonic Tunnel 
is a closed circuit, single-retum atmospheric tunnel. The 
test section is 14.5 ft high by 21.75 ft wide and was 
configured with closed walls and ceiling for this test. A 
more detailed description of the wind tunnel is given in 
reference 1. The computational results obtained in this 
report are compared with experimental data from two 
separate wind tunnel entries with test numbers T421 and 
T437. (Data is in computer storage at the Subsonic Aero- 
dynamics Branch and is available from the author.) 
Unpublished experimental wing-body data from T421 
were compared with computational data for the wing- 
body high-lift Reference H configuration. Unpublished 
data from T437 were used for comparison to the compu- 
tational data for the wing-body with nacelles and 
horizontal tail and wing-body with nacelles and 
horizontal and vertical tail for the sideslip calculations. 
The wind tunnel data (surface pressures and force and 
moment data) from both tests were obtained at a nominal 
dynamic pressure of 0.60 psi, which produced a test 
Mach number of approximately 0.24 with slight varia- 


tions due to temperature changes and a test Re^ of 
approximately 1.4 x 10 5 per in. Pitch sweeps were per- 
formed at angles of attack from -4° to +30°. Sideslip 
angle sweeps at -30° < P < +30° were performed at 
angles of attack of 8°, 10°, and 12°. 

A photograph of the model mounted in the tunnel is 
shown in figure 1. The model is the 0.06-scale Ref- 
erence H high-lift configuration designed by the Boeing 
Commercial Aircraft Group. The model was mounted on 
a close-coupled tandem post mount, which was com- 
posed of a main post and a pitch strut. The same post was 
used for both T421 and T437. A six-component, strain- 
gauge, internal balance measured the forces and 
moments. A spring-loaded sliding seal was fitted to the 
inside of the model to minimize the opening around the 
post mount. A two-point, floor-mounted pedestal support 
system attached near the model aerodynamic center pro- 
vided yaw rotation about the vertical axis for sideslip 
variation and provided pitch rotation about the front sup- 
port point for angle-of-attack variation. The experimental 
data were not corrected for any strut and support interfer- 
ence effects. Forebody grit was applied around the nose 
and along the lower sides of the forebody; there was no 
grit on the wing. A three-view sketch of the model is pre- 
sented in figure 2. Table 1 presents a list of the model 
components and their description. The wing is projected 
to the centerline (see the dotted lines in fig. 2) to deter- 
mine reference area and root chord. 

Table 1 . Geometric Characteristics of Model 


Wing 

Aspect ratio 

Reference area, in 2 

Reference mean aerodynamic 

chord, in 

Span, in. 

Root chord (at the fuselage 

centerline), in 

Outboard break chord, in. . . . 

Tip chord, in 

Location of 0.50c 

Inboard LE sweep, deg 

Inboard TE sweep, deg 

Mid LE sweep, deg 

Mid TE sweep, deg 

Outboard LE sweep, deg 
Outboard TE sweep, deg . . . . 


2.367 

3680.64 

61.937 

93.343 

116.054 

28.282 

7.599 

MFS 126.354, MBL 14.486 

+75.9 

-50.2 

+68.5 

- 8.8 

+48.0 

+10.4 


Vertical tail 

Exposed span, in 13.49 

Exposed aspect ratio 0.864 

Exposed root chord 25.392 

Tip chord, in 5.978 


Horizontal tail 

Exposed span 10.847 in/side 

Exposed aspect ratio 1.297 

Exposed root chord, in 26.752 

Tip chord, in 6.702 
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The pressure ports are located on the Reference H 
model as depicted in figure 3. Though the pressure port 
locations for the upper surface are depicted on the right 
wing and the lower port locations are depicted on the left 
wing, the model had both upper and lower ports on the 
left wing. Tlie experimental pressure distributions at the 
seven designated MFS’s and five designated MBL’s are 
used to correlate with the computational results. 

Computational Grids 

The grids used in this study were provided by the 
Langley-sponsored Geometry Laboratory, which is an 
on-site contracting service. Each grid was generated in 
several stages. First the Initial Graphics Exchange Speci- 
fication data are read into the Integrated Computer-Aided 
Engineering and Manufacturing software system where 
the grid block topology is set up and the coarse Euler 
volume grids are generated. The block faces are refined 
by GRIDGEN2D (an interactive software program used 
to generate grids) to ensure smoothness and orthogonal- 
ity. Additional grid points are distributed in the boundary 
layer for the viscous spacing requirement and the volume 
grids are generated accordingly. 

Several grids of the high-lift Reference H configura- 
tion were generated, all of which had leading-edge flaps 
deflected at 30° normal to the hinge line, 8 l E = 30°, and 
the trailing-edge flaps deflected normal to the hingeline 
at 10°, 5 x£ = 10°. Two of the configurations modeled the 
semispan: the wing-body (configl) and the wing-body 
with nacelles and horizontal tail (config2). The third grid 
modeled the full span wing-body with nacelles and 
empennage (config3) for the sideslip calculations. A 
description of the geometries and conditions solved are 
tabulated in table 2, where a V denotes the presence of a 
component. The wing thickness transitions from a blunt 
leading-edge airfoil to a sharp leading-edge airfoil at the 
planform break, which is located at the 24.35 in. MBL. 
The wind tunnel model makes this transition abruptly, 
resulting in a sharp step; however, the numerical model 
uses a linear interpolation over a region of 0.30 in. to 
smooth this sharp transition region and to simplify the 
grid generation. The leading-edge flaps immediately 
inboard and outboard of the crank are blended together 
so there are no spaces between the flaps. The wind tunnel 
model was also blended together in this region with a 
filler material. The numerical grid did not include the 
wind tunnel walls nor did it model the post support, 
which was used to mount the model in the tunnel; all cal- 
culations were done in free air. A far-field view of the 
surface and volume grid illustrating the overall blocking 
structure, as well as the grid topology, can be seen in fig- 
ure 4. The forebody is a C-0 topology and the wing and 
wake region have a C-H topology. A near-field view of 
the surface is shown in figure 5 and reveals the complex 


blocking structure around the trailing-edge flaps and 
nacelles. The blocking strategy on the lower surface 
around the nacelles is shown in more detail in figure 6. 
The grid geometry on the flaps and flap edges is depicted 
in figure 7, and a close-up view (insert) shows the side- 
of-flap region and flow-through region. An innovative 
multiblock strategy was employed at the trailing edge of 
the wing, which allowed one-to-one matching across the 
trailing-edge flap block interfaces. 


Table 2. Configuration Geometries and Conditions 


Geometry 

Wing- 

body 

Nacelles 

Horizontal 

tail 

Vertical 

tail 

Condition 

configl 

V 




a = 6,8, 10,12,15° 

config2 

V 

V 

V 


a = 8,10,12° 

config3 

V 

V 

V 

V 

a = 8° 

P = 0,6,12,18° 


The full-span surface grid configuration of the Ref- 
erence H 5 le = 30°, =10° with nacelles and empen- 

nage is depicted in figure 8. This grid was generated by 
first coarsening the grid of config2. The reduction in 
points was necessary to ensure that the full-span grid 
could be run in the maximum memory category; compro- 
mises are discussed in the section “Results.” This point 
reduction affected the forebody, the wing upstream of the 
trailing-edge flap hingelines, and the wake regions. All 
other blocks remained the same. The streamwise cuts on 
the wing were decreased by 8, and the grid points in the 
normal directions were reduced by 8. The forebody grid 
was reduced by 48 points circumferentially and 8 points 
in the normal direction. The upper and lower wake 
regions were decreased by 69 points circumferentially 
and 8 points in the normal direction. These adjustments 
resulted in the reduction of over a million grid points for 
the semispan grid. The coarsened grid was then mirrored 
to the other side and point-to-point connections were 
then specified in the longitudinal direction at the fuselage 
centerline. The vertical tail was also incorporated onto 
the geometry. Table 3 summarizes the number of blocks, 
number of grid points, and number of patched and point- 
to-point continuous interfaces for each of the above men- 
tioned configurations. 


Table 3. Number of Configuration Blocks, Grid Points, 
and Interfaces 


Geometry 

Blocks 

Grid points 

Continuous interfaces 

Patched 

Point-to-point 

configl 

14 

3,988,514 

2 

23 

config2 

37 

4,595,343 

38 

83 

config3 

74 

7,085,708 

96 

163 
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Computational Method 

The subsonic flow field about configl, config2, and 
config3 was simulated using CFL3D, This code solves 
Reynolds averaged, thin-layer, three-dimensional, time- 
dependent, compressible Navier-Stokes equations in 
conservation law form. More details about the computa- 
tional method and code can be found in reference 2. The 
interpolated coefficient values between the patched inter- 
faces were determined using RONNIE (ref. 3) and 
PATCHER (ref. 4). The RONNIE code was used for all 
calculations involving configl and config2, and the 
PATCHER code was used for the full-span calculations 
on config3. 

Viscous cell spacing normal to the forebody, wing- 
fuselage-afterbody and wake region was 2 x 1CT 6 c, 
which was required for y + = 1.0 to accurately simulate 
turbulence. A single exponential function was used for 
grid expansion. To minimize the total number of grid 
points, viscous spacing was not specified around the 
nacelle surfaces or at the sides of the trailing-edge flaps. 
The coarse grid spacing at the side of trailing-edge flaps 
is 7.3 x 10“ 4 c. The far-field boundary extends approxi- 
mately 5.8c upstream of the fuselage, 5c circumferen- 
tially away from the center line, and 3.7 c downstream of 
the afterbody. Although the model was strut mounted in 
the wind tunnel, the computations did not simulate the 
strut; all computations done were free-air calculations. 
As previously noted, the experimental data were not cor- 
rected for any strut and support interference effects. The 
conditions of no-slip and impermeability with zero- 
normal-gradient for pressure and temperature are 
imposed at the solid boundaries. Symmetry conditions 
are used along the vertical plane of symmetry for the 
semispan grids. One-dimensional characteristic inflow/ 
outflow boundary conditions are used in the far field. 
Flow-through boundary conditions are specified between 
the trailing-edge flaps and at the nacelle inlet and exit! 

The location of flow transition was not experimen- 
tally determined and therefore unknown. The flow was 
assumed to be fully turbulent for all calculations. The 
field equation turbulence model of Sparlart-Allmaras 
was attempted but the model was unable to handle the 
singular lines in the triangular-shaped regions around the 
trailing-edge flaps. The algebraic turbulence model 
developed by Baldwin and Lomax was therefore used in 
this investigation. For separated flow regions, the method 
of Degani and Schiff was used to determine the proper 
length scale for the turbulence model. The turbulent vis- 
cous terms were calculated only in the direction normal 
to the wall. 

All computations were performed on the numerical 
aerodynamic simulation Cray-2 computer, located at 
Ames Research Center, Convergence was accelerated 


using mesh sequencing and three levels of multigrid. The 
residual and C L history plots are shown in figure 9 for the 
Reference H config2 at a = 10°. The mesh sequencing 
process consists of approximately 1700 iterations com- 
puted on the coarse level, 1500 iterations computed on 
the medium level, and 1000 iterations computed on the 
fine level. An almost 3.5 order of magnitude reduction in 
the residual is observed for the config2 case. Configl 
results showed similar convergence histories. The resid- 
ual and Ci history plots for the config3 case at (3 = 6° and 
a = 8° are shown in figure 10, and a similar mesh 
sequencing strategy was followed. A 3.0 order of magni- 
tude decrease is noted for the config3 case. The algo- 
rithm takes approximately 10 jxsec per iteration per grid 
point. At approximately 37 words per grid point, the 
config3 case requires 262 MW of memory and the 
configl case requires about 170 MW. A typical config2 
case took about 18 hours to converge and config3 took 
about 25 hours. 

Results 

The computational results for the three configura- 
tions are presented below. The general format is presen- 
tation, discussion, and comparison of the computational 
data and experimental data. Comparison of forces and 
moments, as well as surface pressure distributions, are 
discussed for each case. Detailed computational off- 
surface and surface flow visualization are presented to 
help explain the experimental data. 

Configl Results 

The force and moment comparisons for configl at 
6° < a < 15° are shown in figure 11. The computational 
results compare well with experiment, and the trends are 
predicted well. Computational results show that at about 
a = 12°, a vortex-dominated flow develops. The vortex- 
dominated flow causes the nonlinear C L characteristics, 
increase in drag, and nose down pitching moment due to 
greater lift aft of the moment center. A notable difference 
between the computational and experimental results 
occurs at the higher angles of attack (a > 11°) and at 
a = 6°. The computational lift and drag values are higher 
than the experiment; this difference is most likely due to 
the greater influence of the predicted vortical flow. 
Unlike the experimental C m curve, which is relatively 
flat for a < 8°, the computations show a decrease in C m 
at a = 6°. 

Computational pressure distributions are compared 
with experiment in figures 12-21. As discussed previ- 
ously, C p data for the seven MFS and five MBL (fig. 3) 
are plotted. The computational C p distributions on the 
wing and leading- and trailing-edge flap surfaces match 
well with experiment for the lower angles of attack. The 
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comparison between computation and experiment begins 
to show poor agreement somewhere around a = 12°. In 
particular, the outboard pressures on the wing, including 
the crank, show poor comparisons at a = 12° and 15°. 
The surface flow computations indicate a rather massive 
region of separated flow, and the effects of this separa- 
tion are seen as suction peak reductions (see figs. 18 
and 20 at x > 126.86 in. and fig. 21 at y > 28.90 in.). 
These reductions in the suction peaks occur directly as a 
result of vortex dissipation over the wing near the trailing 
edge, and will be presented and discussed later. These 
phenomena do not occur experimentally. 

To show the flow physics on this high-lift high- 
speed civil transport, six crossflow planes showing total 
pressure p t contours are plotted on the Reference H upper 
surface and forebody for a = 6°, 8°, 10°, 12°, and 15° 
(figs. 22-26). The planes shown represent actual grid 
stations and correspond to the following streamwise 
dimensions, which are referenced at the wing leading 
edge: x = 70.2 in., 84.2 in., 103.8 in., 120.2 in., 131.8 in., 
and 154.7 in. At a = 6° and 8°, the flow over the wing is 
attached, except for a small inboard vortex and a tip vor- 
tex (labeled as A and B in figs. 22 and 23). The inboard 
vortex emanates from the wing apex and travels down- 
stream near the wing-fuselage junction. The contour 
lines show a possible vortex forming on the fuselage, but 
ex aminati on of the velocity vectors in that region only 
reveals that the flow is just turning sharply. p t contours 
plotted in crossflow planes at a = 10° are shown in 
figure 24. A more interesting and complex flow field 
develops at this angle of attack. A vortex (C) emanates 
from the forebody and forms over the fuselage. In addi- 
tion to the apex vortex which appears in the two previous 
angle-of-attack cases, a weak elongated vortex (D) sepa- 
rates from the leading-edge flap hingeline and forms just 
outboard of the apex vortex. This vortex first appears at 
about x = 120.2 in. and continues to the trailing edge. 
Another small flat vortex (E) forms at the crank and 
travels along the leading-edge flap until it reaches the 
wing tip and merges with the tip vortex (B). 

The flow field characteristics just described are simi- 
lar for the a = 12° case; however, dissipation of the com- 
bined crank and leading-edge vortices (E) over the 
outboard section of the wing causes extensive flow sepa- 
ration (F and G). At a = 15° (fig. 26), a vortex (A) ema- 
nates from the wing apex similar to the lower angles of 
attack; however, the vortex is larger in shape and magni- 
tude and its footprint traces a more outboard path on the 
wing due to the increased spanwise flow. Just adjacent to 
and outboard of the apex vortex, a leading-edge flap 
hinge line vortex develops (D), which is also stronger 
and more defined than the one predicted at a = 10° 
and 12°. A vortex (C) forms on the forebody and remains 
confined to the upper fuselage, and is similar in shape to 


the one at the lower angles of attack. Just upstream from 
the crank, another vortex (not shown) separates from the 
leading edge of the wing and remains on the flap. This 
vortex does not remain distinct very long since it merges 
quickly with the crank vortex. This combined outboard 
crank and leading-edge vortex elongates and moves onto 
the upper surface of the wing as it propagates down- 
stream and forms vortex H. 

Config2 Results 

The nacelles and horizontal tail were added to 
configl to create config2. The force and moment data are 
plotted with experiment in figure 27 at a = 8°, 10°, 
and 12°, and very good comparisons are seen. The exper- 
imental Q,, C D , and C m curves are very similar to those 
for configl shown in figure 11. It should be noted that 
the wind tunnel model used for config2 had the empen- 
nage attached; however, the computational model only 
had the horizontal tail. The influence of the vertical tail 
on the longitudinal performance and wing surface 
pressures are expected to be minimal so numerical mod- 
eling of the vertical tail was not included. The numerical 
C p distributions are compared to experiment in fig- 
ures 28-33 for a = 8°, 10°, and 12°. A cross section of 
the geometry that shows the location of nacelles and 
trailing-edge flaps are included for the last three fuselage 
stations. This cross section facilitates the interpretation 
of the pressure plots that have the spikes around the 
nacelles and flap edges. Comparison with experiment is 
good for a = 8° and 10°, and in general the C p distribu- 
tions are similar to the Reference H wing-body configu- 
ration, which suggests a similar flow pattern for each. 

The effect of the nacelles influences the surface pres- 
sures close to the trailing edge of the wing since the 
trailing-edge flaps are deflected. As expected, higher C p 
values are seen on the lower surface in the vicinity of the 
nacelles (particularly at x = 143.39 in.). The computa- 
tional Cp results for configl and config2 at x> 135.13 in. 
(figs. 18 and 32) and for y = 28.90 in. (figs. 19 and 33) 
show nonsmooth pressure distributions and poor compar- 
isons with experiment. These nonsmooth computational 
pressures may be due to some flow unsteadiness that 
occurred as a consequence of the flow separation on the 
outboard portion of the wing. The associated computa- 
tional convergence oscillations are more apparent for 
config2 and may be due to grid definition. To accommo- 
date the extra grid points required by the nacelles and 
horizontal tail, fewer points were used in the wake region 
than for the wing-body alone (configl). Hie effect of 
using fewer points in the wake may have influenced the 
upstream wing pressures to a certain extent. 

Off-surface C p contours are depicted in fig- 
ures 34-36 for a = 8°, 10°, and 12°, respectively. 
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Figures 34(a), 35(a), and 36(a) depict results of configl 
and figures 34(b), 35(b), and 36(b) show results for 
config2. Each cross-sectional plane was obtained by slic- 
ing through the volume grid and solution at a specified 
constant x station. These planes are located at x = 1 50 in. 
(cuts across the trailing-edge flaps near the wing tip), 
x= 180 in. (cuts through the fuselage between the wing 
trailing edge and the empennage), and 208 in. (cuts 
through the tip chord of the horizontal tail). 

Analysis of the off-surface flow of configl and 
config2 reveal similar vortical flow patterns at y < 23 in. 
In this region, the nacelles generally reduce the negative 
pressure values in the vortical cores at x = 180 in. and 
affect the locations of the vortices at x = 208 in. Outboard 
of that region, the vortical patterns generated by each 
configuration become more dissimilar. Some of this 
dissimilarity may be due to the difference in the wake 
grid density as previously mentioned. The biggest dis- 
crepancy in the vortical flow patterns between the nacelle 
on and off case is evident for a = 12° at y > 16 in. and at 
x = 150 in. and 208 in., as shown in figures 36(a) and 
36(b). It appears that the horizontal tail has little affect on 
the flow field at all angles of attack when compared to 
the tail-off contours. In fact, the local flow characteristics 
appear to remain almost the same with increasing angle 
of attack. 

Config3 Results 

The results for config3 are presented in fig- 
ures 37-52. As previously stated, config3 modeled the 
full-span high-lift Reference H geometry including 
nacelles and vertical and horizontal tails. Computational 
results were obtained for p = 0°, 6°, 12°, and 18°, at 
a = 8°. Force and moment comparisons between experi- 
ment and computation are shown in figure 37, and fairly 
good comparisons are seen. The experimental values 
remain almost constant from -3° < P < 3° and then a 
gradual decrease in lift occurs over the remaining range 
of sideslip angles. The computational results predict the 
same trend but a more gradual decrease in lift is pre- 
dicted. The numerical Cp, and C m results also compare 
well with experiment; however, since there is an almost 
constant difference between the computational and 
experimental pitching moment, a C^ 0 shift is suspected. 

There are no significant changes in the experimental 
Cj) and C m curves until the sideslip angle exceeds 15°. At 
p > 15°, an increase in drag occurs, which peaks at 
p = 21°. The pitching-moment curve shows decreasing 
slope and transitions to a steep positive slope at 
15° < P < 21°. This change in slope is most likely due to 
an increase in the vortical flow forward of the moment 
reference center (0.05c), which causes a nose-up pitching 
moment and an increase in drag. The measured side force 


coefficient Cy values for p < 15° are linear for that range; 
however, the predicted side forces, which are generally 
underpredicted, depict a nonlinear trend. The numerical 
rolling moment Q and yawing moment C„ compare well 
with experiment. Consistent with the C D and C m curve 
characteristics in the interval, 15° < P < 21°, the Q and 
C„ curves show an increase in slope. A negative Q slope 
is seen for 0° < P < 21° and means that the aircraft is 
directionally stable in roll due to sideslip. Likewise the 
positive C„ curve slope over the same sideslip interval 
indicates that the aircraft is directionally stable in yaw. 
The fuselage and nacelles may contribute to or detract 
from the airplane lateral stability, depending on whether 
the side force of these components acts below the center 
of gravity. The side force acting on the vertical tail 
diminishes the bank angle at the lower yaw angles and 
produces stabilizing characteristics. However, at the 
larger sideslip angles such as those greater than or equal 
to 18°, the flow on the left side of vertical tail stalls and 
produces destabilizing characteristics. 

The C p comparisons for config3 are presented in 
figures 38-47. The spanwise distance plotted on the 
abscissa ranges from -1.0 to +1.0, in which -1.0 to 0 
represents the left side of the wing and 0 to +1 .0 repre- 
sents the right side. Good comparisons are seen for the 
P = 0° case shown in figures 38 and 39. As previously 
mentioned, the grid for config3 was coarser in some 
areas compared to the grid for config2. A detailed 
description of the grids for each configuration was dis- 
cussed in the section “Computational Grids.” A grid 
study was not attempted here since any grid refinement 
would exceed the maximum memory allocation. Since 
the calculations had already been done for config2 at 
a = 8° and p = 0°, the config2 results are compared to the 
config3 results, and the effect of the mesh on the solution 
and flow physics is addressed. The effect of grid density 
on the solution is studied by plotting the C p distributions 
for both cases together as shown in figures 40 and 41. 
The numerical C p distributions are the same though a 
small y shift of the inboard suction peak is predicted for 
config2 at x = 135.13 in. and x = 143.39 in. Most of the 
notable differences occur on the outboard portion of the 
wing near the trailing edge and wing tip. Examination of 
the pressures at the buttline stations in figure 41 shows 
differences occurring at the trailing edge (y = 28.90 in. 
and y = 34.52 in.) as well as near the wing tip at 
y = 42.0 in. 

The C p comparisons for p = 6°, 12°, and 18° for 
config3 are shown in figures 42-47. Note the asymmetric 
C p distributions characterized by higher suction peaks on 
the right wing compared to the left wing, which causes a 
negative rolling moment (see fig. 37). The increased 
effect of the fuselage, nacelles, and trailing-edge flaps is 
evident in the pressure plots as P increases through 12° 
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and contributes to increased lateral stability (steeper Q 
and C„ curve slopes are evident in the interval 
15°<P<21°). At p = 18°, these components have an 
even greater effect on the surface flow, and very 
erratic and unsmooth pressure distributions are seen at 
y/(b/2)i = 143.39 in. and 150.00 in. The effect of the 
fuselage on the pressure distributions is seen in figure 46, 
which shows a dramatic increase in negative C p around 
y/W)i = 0. 

While the pressure comparisons are generally good 
for p = 6°, 12°, and 18°, some of the predicted suction 
peaks deviate in location or just are not seen in the 
experimental results. For example, the predicted suction 
peak at y/(b/2)i = -0.5 (x = 135.13 in. and x = 143.39 in. 
shown in fig. 42) appears to be inboard of the experimen- 
tal peak, which is located at y/(b/2)i ~ -0.65. Another 
predicted suction peak, which occurs at x = 143.39 in. 
and y/(b/2)i = -0.70 at (3 = 12° (fig. 44), is not verified by 
the experimental data at all. At P = 18° the numerical 
results in the region of the wing leading edge did not 
match well with experiment as shown in the first three 
MFS’s in figure 46. Comparisons at the MBL stations 
were good, and it is noted that the suction peak that con- 
sistently occurred at x/c; = 0.65 at buttline station 42.00 
for the entire range of sideslip angles was not computa- 
tionally predicted. 

Additional details of the flow around config3 at vari- 
ous sideslip angles are seen in figures 48-52. All cases 
were run at a = 8°. For each sideslip angle, the flow char- 
acteristics are first examined by plotting pressure con- 
tours on the upper surface as shown in figures 48-50. 
The flow about the empennage is extracted from the 
numerical solution at x = 200 in., which bisects the verti- 
cal tail at the tip and cuts through the horizontal tail just 
upstream of its tip. The off-surface C p contours and 
velocity vectors are shown in figures 48 and 51, respec- 
tively. To further interpret this empennage flow, stream- 
lines on the left side of the vertical tail and upper aft 
portion of the fuselage are plotted in figure 52. 

Several noteworthy flow characteristics develop as 
the sideslip angle is increased. The suction pressures on 
the right side of the wing increase, which is seen in the 
surface pressure contour plots (figs. 48-50). Examination 
of the off-surface pressure contours and velocity vectors 
about the empennage show the increasing effect of the 
crossflow with increasing sideslip angle. The pressure 
contours (fig. 51) plotted on the left side of vertical tail 
for P = 1 8° reveal a very disorganized pattern of vortical 
flow, and the velocity vectors (fig. 52) show increased 
circulation in this area, as well as the formation of 
additional vortices. This increased circulation contributes 
to the extensive region of stalled flow. The surface 
streamlines, shown in figure 53, reveal a highly three- 


dimensional flow field characterized by complex vortical 
flows and extensive separated regions. For all sideslip 
angles, a primary vortex separates from the sharp leading 
edge of the vertical tail and reattaches near the trailing 
edge of the tail. At P = 12°, the flow on vertical tail 
reveals two saddle points, which develop as a result of 
the surrounding flow converging to one node or point. 
Secondary and tertiary separation lines are also evident 
from this figure. At p = 18°, the flow on the left side of 
the vertical tail reveals a rather extensive region of 
stalled flow. The streamlines confirm the existence of the 
stalled flow and a saddle point (larger than the ones seen 
at P = 12°) develops. As the sideslip angle is increased to 
21° and beyond, die increased effect of the stalled flow 
on the tail may contribute to the destabilizing characteris- 
tics, which were discussed previously (fig. 37). This type 
of separated flow appears to contribute to the increase in 
drag and pitching moment that occurs for the interval of 
15°<P<21°. 

Conclusions 

The subsonic flow field about several high-lift 
(leading-edge flap deflection angle of 30°, trailing-edge 
flap deflection angle of 10°) Reference H configurations 
was successfully predicted using CFL3D, a Reynolds 
averaged, thin-layer Navier-Stokes flow solver for struc- 
tured grids. An innovative multiblock strategy was 
employed at the trailing edge of the wing, which allowed 
one-to-one matching across the trailing-edge flap block 
interfaces. Three versions of the Reference H model 
were simulated, each possessing different components 
attached to the basic wing-body geometry. The three 
cases analyzed were the wing-body, wing-body with 
nacelles and horizontal tail, and wing-body with nacelles 
and horizontal and vertical tails. The first two configura- 
tions modeled the span to the symmetry plane, and 
solutions were obtained at angle of attack 6° < a < 15° 
and zero sideslip. The third configuration modeled the 
full-span geometry in order to analyze the flow physics 
and aerodynamic performance at different sideslip angles 
of 0° < P < 18° at angle of attack 8°. All calculations 
were performed in free air at wind tunnel conditions of 
Mach number 0.24 and Reynolds number 1.4 x 10 5 per 
in. The post-mount model support system used for the 
test was not computationally simulated. The applicability 
of the structured Navier-Stokes method is demonstrated 
by comparing the predicted results with the experimental 
surface pressures and force and moment data. 

The longitudinal performance characteristics were 
well predicted and very good force and moment compar- 
isons are seen. The pressure distributions show generally 
good comparisons and trend at the lower angles of attack; 
however, comparisons diverge when the computational 
flow on the outboard portion of the wing becomes 
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extensively separated, which occurs at about 12° angle of 
attack. Improvements in the solution may be obtained by 
limiting the Degani-Schiff search to the boundary layer 
in the Baldwin-Lomax turbulence model as implemented 
in CFL3D. Agreement may also be improved by specify- 
ing the experimentally determined transition location 
from laminar flow to turbulent flow. Not knowing this 
information would affect the computational results since 
transition determines the onset of separation and drag. 

Postprocessing of the numerical solutions provided 
detailed surface and off-surface flow quantities not 
obtained in the wind tunnel. The plotted off-surface pres- 
sure contours of the computational vortical flow field 
were very useful in interpreting the wind tunnel data and 
provided a comprehensive three-dimensional under- 
standing of the flow. As is evident from the plotted off- 
surface flow contours, a very complex multiple vortical 
system develops at the higher angles of attack (a > 12°). 
Also, at the higher angles of attack, the flow on the 
inboard wing turns spanwise onto the outboard panel. 
This spanwise effect further complicates the highly vorti- 
cal flow since it impacts the boundary layer and causes 
extensive separation on the outboard section of the wing. 

The lateral performance characteristics ^howed 
fairly good quantitative comparisons, and the overall pre- 
dicted trends in the yawing moment, rolling moment, and 
pitching moment were well predicted. The computational 
surface pressures that are extracted at constant spanwise 
stations reveal asymmetric pressure distributions charac- 
terized by higher suction peaks on the right wing com- 
pared to the left wing. The side force on the vertical tail 
at the lower yaw angles produces a stabilizing effect. 
However at the higher sideslip angles, an extensive area 
of stalled flow forms on the left side of the vertical tail 


and contributes to the destabilizing characteristics of the 
model. 

The high-lift high speed civil transport work is ongo- 
ing in support of the High Speed Research program, and 
future work will be extended to updated configurations in 
the downselect process. Because of the complexity of the 
multiblock gridding, as well as the total number of grid 
points required to adequately resolve the high-lift flow, 
this method would not be considered for design. How- 
ever, it does well as an analysis tool and future improve- 
ments in the CFL3D code should be made to broaden its 
application to include more accurate simulation of multi- 
ple vortical flows and boundary layer separation. 


NASA Langley Research Center 

Hampton, VA 23681-2199 

March 17, 1998 
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Figure 1. Boeing 0.06-scale Reference H model with deployed leading- and trailing-edge flaps in Langley 14- by 
22-Foot Subsonic Tunnel. 
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Figure 3. Schematic representation of upper and lower surface pressure port locations for 0.06-scale Reference H 
model. All pressure ports were located on left wing (below centerline) on model. 
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Figure 6. Close-up view of multiblock arrangement around nacelles. 
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Figure 13. Concluded. 
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Figure 15. Concluded. 
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Figure 17. Concluded. 
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Figure 19. Concluded. 
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Figure 21. Concluded. 
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Figure 28. Spanwise distributions of experimental and computational C p for Reference H config2 at M = 0.24 
and a = 8°. 
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Figure 33. Concluded. 
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(b) Config2, a = 12°. 
Figure 36. Concluded. 
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Figure 41. Concluded. 
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Figure 45. Concluded. 
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Figure 47. Concluded. 





























Figure 53. Surface streamlines on left side of vertical tail and fuselage for config3 at a = 8' 
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